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^ ■ Abstract 
o 

Q>^ \ We investigate ground-state and thermal properties of a system of non-relativistic bosons 

2 . interacting through repulsive, two-body interactions in a self-consistent gaussian mean-field 

s 

. approximation which consists in writing the variationally determined density operator as 
O ■ the most general gaussian functional of the quantized field operators. Finite temperature 



X 



results are obtained in a grand canonical framework. Contact is made with the results of 



H ' Lee, Yang and Huang in terms of particular truncations of the gaussian approximation. The 
full gaussian approximation supports a free phase or a thermodynamically unstable phase 
when contact forces and a standard renormalization scheme are used. When applied to a 
Hamiltonian with zero range forces interpreted as an effective theory with a high momentum 
cut-off, the full gaussian approximation generates a quasi-particle spectrum having an energy 
gap, in conflict with perturbation theory results. 
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1 Introduction 



The non-relativistic interacting Bose gas is certainly among the most extensively investigated 
problems of many-body physics. Early interest in this problem has been strongly motivated 
by the low temperature properties of helium |l|]. It has been sustained, on the theoretical 
side, by a persisting elusiveness of the deeper nature of the condensation process and, 
on the experimental side, by developments which led to therecent observation of the more 
artificial Bose condensates 0-0. 

Much of the persisting uneasiness about the condensation process of a nonideal bose gas 
hinges on the fact that the body of available results relies mainly on perturbative methods, 
a circumstance which tends to be considered as a liability. Non perturbative approximation 
schemes have been developed, on the other hand, in connection with the problem of self- 
interacting, relativistic Bose fields which became relevant e.g. for inflationary models of 
the universe This problem has also been used as a testing ground for non perturbative 
methods which could then be applied to more complicated systems of interacting fields. In 
this context, the gaussian variational approach |^ has received considerable attention, also in 
view of its relation to extended mean field methods traditionally employed in non relativistic 



many-body physics M, 0, 10, 11 



The purpose of this paper is to describe an application of the gaussian variational ap- 
proximation to the much studied problem of interacting non relativistic bosons in order to 
bring about a direct confrontation of this approach with at least part of the available results 
for this problem. We use a formulation |p which is very close to the standard language 
of non relativistic, extended mean field approaches of the Hartree-Fock-Bogolyubov type, 
which can on the other hand be directly related to the methods adopted in connection with 
the Schrodinger representation of quantum field theory |10|. Although for simplicity we 



restrict ourselves to the case of uniform systems, extending the formulation to finite, in- 
homogeneous systems such as those actually realized in the recent alkeli atom experiments 
is completely straightforward ||^ using well known many-body techniques [|12[. The inclu- 
sion of finite temperature effects is also straightforward in the formulation we use, so that 
thermodynamic properties can be studied rather easily. We show that several of the old 
results can be retrieved from truncated versions of the full gaussian approximation, including 
ground state energy and phonon spectrum. This last feature is lost when one adheres to the 
full gaussian approximation, notwithstanding the fact that it is supported theoretically by 



the Hugenholtz-Pines theorem |[T^, |1^ to all perturbative orders, a feature which has been 
noted long ago by Girardeau and Arnowitz ||15|. The phonon spectrum can be recoverd in a 
generalized RPA treatment or by using the functional derivative method developed by 
Hohenberg and Martin Jl^ and recently re-discussed by Griffin [T^. Problems of thermo- 
dynamic instability develop for contact forces when one uses a renormalization procedure 
in which the coupling constant is made to approach zero from it negative values, which has 



been considered in a relativistic context e.g. by Bardeen and Moshe |jT9[ and which has also 
been extensively discussed by Stevenson as the "precarious" renormalized theory |^ 



and more recently by Tiirkoz and by Kerman and Lin |TT 



The formulation which we adopt is reviewed in Section 2, where the formal, finite tem- 
perature equilibrium solutions for a contact repulsive interaction are obtained. A truncated 
version of the gaussian variational equations is analyzed in Section 3. In Section 4, we deal 
with the full gaussian approximetion. Two different prescriptions for dealing with the di- 
vergences are considered. The first one involves a renormalization scheme related to that 
proposed by Stevenson. As a second, alternative scheme we treat the Hamiltonian with con- 
tact interactions as an effective theory to be implemented with a fixed cut-off in momentum 
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space. Numerical results are given in this case for properties of the different phases and 
phase equilibria. Section 5 contains our conclusions. 

2 Thermal gaussian approximation 

We consider an extended, uniform and isotropic system of non-relativistic, interacting, spin- 
less bosons described by the Hamiltonian (in momentum representation with periodic bound- 
ary conditions in volume V) 



H = e(k)a^y.aT: -\ — — a]^ ji]^ -Hcac il] 

^ ^ ' k k ' 2Y ^ fci+q fc2-g Ki ^ > 



2V 

k k\k2q 



where e(/c) = is the free particle kinetic energy, and the contact repulsive (A > 0) 
interaction between a pair of particles is \5{f — r'). The field operators satisfy standard 
boson commutation relations. Within a grand canonical description the state of the system 
is described by the density operator 

Z 

where Ti = H — fiN, N being the boson number operator, and Z the grand canonical 
partition function 

Z = Tre-^^. 



Following the variational approach of Balian and Veneroni |^T|, we look for extrema of the 
object 



/(M) = Tre 



-M 



under the constraint M — (31-i = 0, which is taken into account through the introduction of 
a Lagrange multiplier matrix B. This leads to the variational problem 



5$(M, B) = 5Tr e"^^ + B{M - PU) 



0. 



(2) 



Variation of M gives B = e and elimination of B from Eq. |^ leads to the M-dependent 
object 



$(M,5) ^(M) = Tr e"'"(l + M - 



(3) 



Balian and Veneroni show that Z > \E'(M) for any M. Defining z = Tre we can write 



where now TrjFg = 1. Substituting this in Eq.|^ and varying z we find the variational 
expression for the grand potential Q 



Q = -^In Z <Tr[{n + KT In J^o)J^o)] 



(4) 



where J-'o is an arbitrary density with unit trace, and we used the notation 1/(3 = KT. In 
Eq. I we can in particular identify an entropy factor as 5*0 = —K Tr [JFq In jFg]. 

The most general gaussian approximation consists in adopting for M an ansatz of the 



form [22 



M ^ Y.[^rk'<^hk' + iBru.a\al + h.c.) + (Qat + h.c.)] 

kk' 

where the matrix A^:^, is hermitean. The quadratic form appearing in M can be diagonalized 
by a general canonical transformation of the Bogolyubov type, which amounts to changing 
to the natural orbital representation of the extended one boson density corresponding to the 
gaussian density JFq [||, |^. The uniformity and isotropy assumptions we make allow us to 
restrict this general ansatz so that is diagonal, -B^p vanishes unless k = —k' and both 
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of these matrices and the depend only on the magnitudes of the momentum vectors. The 
diagonahzation of the quadratic form is achieved in this case by defining transformed boson 
operators as 

m = ^Ih + y*Ak 

where 

and we have used the isotropy of the uniform system to make the c-number transformation 

— * 

parameters Xk-, Vk and dependent only on the magnitude of A;. In order for this trans- 
formation to be canonical we have still to impose on the Xk and yk the usual normalization 
condition 

\xk\^ - \yk? = (5) 
The trace-normalized gaussian density operator is now written explicitly as 




Straightforward calculation shows that 

TrivlVk'^o) = ^kSkk' 

so that the i/k are positive quantities corresponding to mean occupation numbers of the 
77-bosons. One also finds that 
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with Ak — xlFk + yiri so that non vanishing values of the correspond to coherent 
condensates of unshifted, Bogolyubov transformed bosons. We again invoke the system 
uniformity to impose 

in the calculations to follow. 

It is important to note that the truncated density in general breaks the global gauge 
symmetry of H which is responsible for the conservation of the number of a-bosons. The 
quadratic dispersion of the number operator N — J2k ^^^^k ^^^^ state can in fact be obtained 
explicitly as 

{N')-{Nf = 2|ron|a;o|% + |2/or(l + ^/o)]-2r;'a;oi/^(l + 2i/o)- 

-2Tlyoxl{l + 2uo) + irol' + Y.[\xk\^i^k + (1 + i^kMl^] + 

k 

+ Y.iW^'klW + (1 + + kferiyfer(i + 2i^ifc)'}- 

k 

Furthermore, mean values of many-boson operators taken with respect to J-'q will contain no 
irreducible many-body parts, so that the replacement of T by amounts to a mean field 
approximation. The states described by JFq have therefore to be interpreted as "intrinsic" 
mean field states. 

An important simplification which occurs in the case of stationary states such as we 
consider here is that the transformation parameters Xk, Vk and Fq can be taken to be real 
and we can use a simple parametric representation that automatically satisfies the canonicity 
condition. It reads 
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Xk = cosh ak, = sinh cTfc. 
It is then straightforward to evaluate the traces involved in Eq.|^ to obtain 



(7) 



^ < 5] e(fc)-/i + 



V 



[1 + 2z/fc) cosh2crfe — 1 



2 , 

2V 



[1 + 2z/fc) cosh2crA: — 1 



-^E(l + 2^^)sinh2a, + A E 

k { k 

+ ^{E(l + 2^^-)sinh2a,f 

k 

-KTY^[{1 + Vk) ln(l + z/fc) - z/fclnz/fc]. 



In a similar way the number fixing condition Tr [jFgA^] = (A^) evaluates to 



{n) = ti+y: 



'1 + 2vk) cosh 2crfc — 1 



(9) 



2.1 Formal equilibrium solutions 

Equations determining the form of the truncated density JFq appropriate for thermal equilib- 
rium are in general derived by requiring that f2, eq.^, is stationary under arbitrary variations 
of Fq, CTfc and z/^. Variation with respect to Vq gives the gap equation 



2A A 

To <! —Tl - 2/i - - X] [(1 + 2z/fc)(sinh2(T, - 2 cosh 2^^) + 2] 

k 



(10) 

which, besides the trivial solution Fq = 0, may also admit a solution with a non vanishing 



value of Fq obtained by requiring that the expression in curly brackets vanishes. This solution 
involves the number constraint, Eq.^ in addition to the values of Vk and ak-, which are 
determined by the remaining variational conditions on f2. In order to simplify the algebraic 



9 



work involved in the study of this class of solutions it is convenient to use the number 
constraint Eq.|^ to eliminate Fq from the right hand side of Eq.^] which then assumes the 
form 

n<F- i2{N) 

with (A^) given by Eq.O. This identifies a free energy F as 



[1 + 2z/fc) cosh 2o"fc — 1 



it 



[1 + 2z/fc) cosh 2o"fc — 1 



+ ^ + 2^^^') sinh2(Tfc' 

k,k' 



+ 3T7{E(l + 2^^0sinh2afe}2 
oV 

k 



[1 + 2vk) cosh 2o"fe — 1 



+ 



Xp^V 



-KTY^[{1 + Vk) ln(l + z/fc) - Vk In z/^] 



Extremizing F by setting derivatives with respect to CTfc and Vk equal to zero one gets 



and 



where 



tanh2afc= (12) 
Ap - 2V + 2^^^) cosh 2(Tfc - 1] - ^ (1 + 2i/fc) sinh 2a, 



e{k) + Ap - 2^ + 2z/fe) cosh 2ak - 1] + 2V (1 + 2z/fc) sinh 2ak 



^k = (13) 

{exp[VA/KT]-l} ^ ' 



A = e{kf + 2\e{k)x 
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^{P- + 2'^fe)cosh2afc - 1] + ^ X] + 2^^^) sinh 2afc} 

k k 

+4A2{p - ^ E[(l + 2^fc) cosh2(T, - 1]} X 

it 

x^E(l + 2'^fc)smh2a, (14) 
As for the trivial solution Fq = 0, we differentiate Q as written in EqM and get 



-^Efc(l + 2t/fc)sinh2afc 
e(fc) - /i + ^ + 2z/fc) cosh 2cTfc - 1] ' 



tanh2afe = ^^aT. r/T . 7T- (15) 



Finally, we stress that the results obtained in the present section are largely formal, as 
they involve divergent sums. In order to allow for the derivation of thermodynamic properties 
of the different phases they must therefore be supplemented by suitable regularization and 
renormalization procedures. These will be discussed in sections 3 and 4 below. 

3 Independent 77-bosons and dilute 
system limit 

In this section we consider two different truncation schemes of the gaussian variational equa- 
tions which lead to well known results for the low temperature properties of the interacting 
boson system. The first and most drastic truncation of the gaussian approximation consists 
in neglecting all terms representing interactions between 77-bosons. The result is entirely 
trivial in the case Fq = since this implies cxfc = 0. All effects of the interaction are thus 
discarded, giving us just the ideal gas results 



{exp[<p^] - 1} ^^^^ 
where fi is determined by the number constraint 
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P-^^k {exp[^]-l}- ^ ^ 

As for the solution corresponding to Fq 7^ 0, discarding interaction between r^-bosons 

amounts to dropping all double sums in Eq. (9). The variational conditions on 0"^ and v^. 

appear then as 



tanh2(7fc^ '^^ (18) 

e[k) + Xp 

and 



= , , = . (19) 

When using these results for calculating F we replace sums by integrals and introduce a 
cut-off A in the range of integration over momenta. Ignoring terms that vanish in the limit 
A — > 00 we get 



F 
V 



Xp" X^mkp" %rrv'l''X''/''p^''' 



2 At^'^Iv' 

KT rA 



+ 



+ 



27r2 



j A;^ In < 1 — exp 



e{k)^ + 2e{k)Xp 
KT 



dk 



(20) 



which shows that we get a linearly divergent term proportional to A^. This term is a leftover 
of the non-normal ordered terms of H which involve two r^-operators. It is therefore a 
direct consequence of the non-trivial nature of the Bogolyubov transformation, and can be 
compensated by introducing the Fermi pseudo-potential, where the contact interaction is 
replaced by 
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Vp{r-fl) = A 



d 



|r — r/\)6{f — r/)] 



(9(|r — rf)\ 
= \6{r-ff) + {\r-ff\)\6'{r-ff) 

which can be related, in the case of dilute, cold systems to the two boson scattering length 
a through 



A 



m 



The regularized ground state energy becomes 



V 2 I 15 7rV2 



and the chemical potential appears as 



Ap + 
1 



4^3/2^3/2;^5/2 
3rf 



27r2 Jo 



Xe{k)edk 



exp 



KT 



^e{ky + 2e(A;)Ap 



(21) 



These are just the results obtained by Lee, Huang and Yang ||23|, by Beliaev ||2^ and by 



Hugenholtz and Pines under the assumption of a macroscopic (c-number) occupation of 
the zero momentum mode. Note however that in the present formulation this is replaced by 
the coherent condensate associated with Fq. In addition to this we still have non-coherent 
occupation of the zero momentum mode as given by the limit of z/^, Eq. (17), for k ^ 0. 
For small, non-zero temperatures this diverges as 1/k and therefore does not contribute to 
the density of the system. Finally, the values obtained for the chemical potential in the 
condensate (Fq 7^ 0) and non-condensate (Fq = 0) phases indicate the instability of the 
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former. This may be seen very easily from the fact that /i is always positive in Eq. (19) and 
always negative or zero in Eq. (15). 

A possible way to circumvent this drawback in the framework of the gaussian approx- 
imation is to perform a somewhat less drastic if more delicate truncation of the complete 
variational expressions, which relies on the fact that we are working with dilute systems. 
There is little change in the case of the Fq 7^ phase. Eqs. (7) and (16) give the depletion 

8(aV)'/' 

In the hmit of a dilute system one has (a^po)^^"^ ~ (a^p)^/^ <S 1 so that one may replace p by 
Po in Eq. (16). With this replacement the excitation spectrum of the truncated Hamiltonian 
is related to the density of the condensate as 

w{k) = -y^(^j2T2Ap^^ 

As for the phase with Fq = 0, the occupation numbers Uk are obtained by keeping 
(jfc = but with no truncation. The chemical potential p is determined using the constraint 
condition (7) that in this case reads 

1 r'^ k^dk 

Since the truncation has been avoided in the case when Fq = 0, this phase is no longer treated 
as a free bose gas, and calculation shows that at T = the condensate phase now appears as 
the stable one, since its chemical potential p — Xp is lower than that of the non condensed 
phase, for which p — 2\p. It is apparent that in this derivation the two phases are not 
treated on the same footing as far as the truncation is concerned. However the truncation of 
the two body terms has a quite different significance in each of the two cases. The rationale 
for this procedure rests in fact on the expectation that the most relevant effects of the two 
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body interaction are incorporated to the result via the symmetry breaking processes which 
take place in the treatment of the Fq 7^ phase but not in that of the Fq = phase. 

4 Handling the full gaussian approximation 

In this section we examine two different prescriptions for dealing with the divergences of 
the complete gaussian approximation developed in section 2. The first prescription involves 
a renormalization scheme similar to the one proposed by Stevenson [^] in the context of 
the relativistic 0^ theory under the name of "precarious theory" , in which the bare coupling 
constant is made to approach zero by negative values. It represents an attempt at sticking 
as much as possible to the dynamics of contact interactions as such in the present context. 
Although successful in removing the divergences in a consistent way, this scheme will be 
shown to lead to a thermodynamically unstable system. We therefore consider also, as a 
second prescription, the simple alternate "effective theory" scheme in which a fixed cut-off 



is introduced in momentum space in the spirit of the work of Amelino-Camelia and Pi 
4.1 Contact forces: precarious theory 

In order to reduce unessential complications to a minimum we restrict ourselves in the 
following development to the properties of the system at T = in the phase corresponding 
to Fo 7^ 0, since the extension to T 7^ involves no additional divergences. The gap equation, 
Eq. (8), combined with the number constraint, Eq. (7), gives 



/X = Ap+ ^^[cosh2afe - 1] - A^sinh2(Tfc. (22) 

k k 



and from Eqs. (10) to (12) we get 
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tanh 2a, = ^ ' ^ ^/t^"^?^ ' ' f ^/ ""^ '"^^ . (23) 
e{k) + Ap - 2y Efc[cosh2(Tfc - 1] + 2y EfcSinh2crfc 

This is an implicit equation for the transformation parameters a, which appear in the sums of 
Eqs. (20) and cut-off A and, again neglecting contributions that vanish in the limit A oo, 
make the ansatze 



^X^sinh2^Tfc = a + /3A (24) 

k 

and 



^^(cosh2afc-l)=7 (25) 

k 

where a, (3 and 7 are assumed to approach finite values in the limit A 00. We next 
introduce a renormalized coupling constant as 



^ (26) 



so that when A — 00 the bare coupling constant A approaches zero from negative values 



n|, |T0|, gg. Eq. (21) becomes 
tanh 2ak 



Arp— A,.7— ArQ— Ar/3A 

_ M+f + 0(A- 



e(fc) + A.P-A.7+WA./3A e(fc) + M + f + 0(A-2) 



with 



M 



m 

.2*2 o_2fc2. 



27i'h\ 2Ti'h'(3^ 

Q = (p-7-a r ), 

m Arm 

R = p-7 + a + ^-^. 27 

m Arm 
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This makes good sense for all k and for any A provided M is positive. We can also obtain 
sinh 2(Tfe and cosh 2(7^ in terms of M, Q and R to evaluate the sums of Eqs. (22) and (23) 
up to terms C(A~^) with the results 



1 /-A k^(M+Q-) 
a + (3A = --r^ ^ =dk 

47r2Jo ^e{ky + 2{M + ^)e{k)+G 

_ m^/^M^/^ mMA 
2^~2^' 

which is consistent with the expression for M in Eq.(25), and 



47r'-^o [^e{ky + 2{M+Q;)e{k) + G j 3rf 
where we defined 

_ 2M(i? - g) i?2 _ g2 

A +^V^- 

It is easy to see that for M = we obtain the ideal Bose gas results (free theory). For 

M > 0, in order to evaluate F (which for T = reduces to the ground state energy Eq ) , we 
also need 

1 ^ ,,,, , ^ mM^A mMQ Am^/'^M^/'^ 
- Ee(^)(cosh2., - 1) = + 

k 

Taking these results to Eq. (9) (with i/^ = for T = 0) we see that the remaining linearly 
divergent terms cancel and we are left with 

^ = — = ^ Mp (28) 

which identifies the chemical potential as = —M. Note that M > now implies < 0. 
Finally, in order to relate ji (or M) to the density of the system we evaluate the grand 
potential, Eq. (6), with the result 
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n(r^o)^?!^^!W^_^ (29) 



and use the relation = —{^)t,v to obtain 



^ 4m3/2(-^ )3/2 

With the condition /i < when A oo, we get a phonon-hke spectrum as in Sec. 3 



p = - + — ^ri73 — • (30) 



E{k) = ^e{kY -2ne{k). 

We can also take the appropriate derivative of F (or, equivalently, evaluate —VL/V using 
Eq. (27)) to get the pressure at T = as 

and it is easy to show, using Eqs. (26), (28), (29) and p > 0, that ^ and EqIV are always 
negative. This shows that the renormahzed theory is thermodynamically unstable. As 
a matter of fact, this instabihty should be seen as the non-relativistic counterpart of the 



"intrinsic instability" pointed quite some time ago by Bardeen and Moshe [IS] in their 
analysis of the phase-structure of the relativistic Xip^ theory in four dimensions, which has 



also been raised recently in the context of multi-field theories [E^. It underlies also the 



"precariousness" of this theory in Stevenson's treatment ||20[| . 

Another possible phase to be studied corresponds to a solution where Fq = 0. In order to 
account for the fact that in this case a macroscopic occupation z/q may develop, we separate 
the A; = contribution in the sums 



^ Y: (cosh 2ak-l)^c + ^Yl (cosh 2a, - 1) (32) 

k k^O 



18 



and 



— — sinh 2a k 



SO that Eq. (13) becomes 



tan 2crjk 



-Ad - ^Efc^oSinh2(7fc 
e{k) - ii + 2\p 



Using the ansatze (22) and (23) we obtain 



tanh 2(7fe 



— Ard— ArQ— Ar/3A 



with 



e{k) - /X + 



2Xrp 



TT 



M+^ + C(A-2) 
e(A;) - /X + f + O (A-2) 



(33) 



(34) 



(35) 



M = 

Q = 

R = 



+ 



m 



m 



m 



-P- 



(36) 



If we examine 



1 



7 = 



47r2 Jo 



e(A;) — /i 



(37) 



Y[e(A;) - /i]2 - M2 

we see that the only acceptable solution is M = /j. If M = it follows that (7^ — 0, leading 
to a possible solution c — d — VQ/V which corresponds essentially a free Bose gas. If, on the 
other hand, M 7^ we are again led to a thermodynamically unstable situation. 
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4.2 Effective theory 

The renormalization prescription of the preceding subsection tells us that if is positive 
and held fixed when A ^ oo we must have A — > 0_ which results in a thermodynamically 
unstable theory that we want to avoid. Since a positive value of A leads to a trivial theory 
with Xr — when the limit A — > oo is taken, we next consider the results that one obtains 
for the non-trivial effective theory in which A is kept fixed at a finite value allowing for both 
Ar > and A > 0. This requires that one must have 

XrinA 

and implies a finite resolution 0{A~^) in configuration space. The restricted momentum 
space also implies that the validity of the results will be restricted to temperatures which 
are low in the scale 

2Km' 

Because in the derivations given in the preceding subsection all terms that vanish in the limit 
A ^ oo were neglected, we use here the general expressions given in Section 2. Eq. (24) will 
be used with a fixed value of A to relate the bare and effective coupling constants A and A^ 
respectively. 

Consider first the condensed phase, Tq ^ 0. The sums appearing in Eqs. (9), (10) and 
(11) are now finite, and we therefore define 

^ = ^X^'(l + 2i^fc)sinh2(7fe 

k 
k 
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k 

where the primed sums are restricted to \k\ < A. Eqs. (8), (10) and (12) appear then as 



fi = Xp-XA + XB (38) 

X[p-B-A] 

tanh2(7r: = \, „ ^ 

e{k) + X\p-B + A] 

and 

A = e{kf + 2Xe{k)[p - B + A] + 4X^[p - B]A. 

This last quantity determines through Eq. (11) and shows explicitly the presence of an 
energy gap in the quasiparticle spectrum which is related to the quasiparticle interaction 
term 2XpoA (we use p — B = po, the density of the coherent condensate). The free energy F 
and the pressure P are given by 

F = CV - Xp{A - B)V + ^{p'^ + A"^ - B'^)V + XABV 

-KT^[{1 + Uk) ln(l + Uk) - Uklnuk] (39) 

k 

and 



XB' Ap2 XA^ 1 ^/ r- 

k 

KT 

-Y^\n{l - expi-VA/XT]} (40) 



V 

k 

After going to the continuum limit, we can calculate A and B for given values of p and 
T by solving numerically the system of equations 
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Finally, in terms of the values of A and B we can calculate C and the thermodynamic 
functions also numerically. 

As for the case Fq = with the sums restricted to the cut-off, it's easy to see that the 
only solution is Ufc = 0. In the continuum limit, 

^ ~ 27r2 Jo exp[(e(A;) - + 2Xp) /KT] - 1 ^ ' 

which serves to determine /i. The pressure is in this case given by 

KT 

P^Xp-"-—- e ln{l - exp[-(e(A;) -p + 2Xp)/KT]}. (42) 

ZTT^ JO 

4.3 Numerical results 
4.3.1 Cut-off dependence 

In order to perform numerical calculations we used a system of units such that the Boltzmann 
constant K = 1 and c = 1, so that energies are expressed in degrees Kelvin. We use 
m — A X 10^^ °K and = 50 °KA^. In Fig. (1), we show the dependence of the ground 
state energy, for a density p = 0.01 A ^ as the cut-off is varied in the interval 1.0 A ^ < A < 
5.0 A ^. For the calculations described below we take A = 3.0 A ^. According to Eq. (24), 
this corresponds to a bare couphng constant A = 119.124''irA^. For densities larger than 
the quoted value one starts having stronger cut-off dependence, so that the present scheme 
is, in this sense, also limited to rather dilute systems. 
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4.3.2 Phase transition 

Results of some sample calculations using the various expressions of the preceding section 
are shown in Figs. (2) to (5). For the range of temperatures and densities considered here, 
the adopted value of the cut-off momentum A essentially saturates the values of the finite 
momentum integrals involved. The density p is plotted in the figures as the number density 
in units of A ^ . 

In order to study the properties of the condensed (po > 0) and of the non-condensed (po = 
0) phases, as well as their coexistence, we show in Fig. 2 the chemical potential ii calculated 
as a function of the density p in each of these two cases, for a constant temperature fixed 
as T = 2°K, using Eqs. (38) and (41) respectively. The condensed phase has two different 
branches with different chemical potentials in the range p >~ .7x 10~^A ^ and ceases to exist 
for lower values of the density. The chemical potential of the non-condensed phase increases 
monotonically approaching the upper branch of the condensed phase asymptotically from 
above. It is interesting to note that these two solutions essentially merge at a density where 
the non-condensed solution develops a macroscopic occupation uq of the zero momentum 
state, and that the solution involving instead a coherent condensate (po > 0) leads to a 
lower value of In order to decide about the thermodynamic stability of the various 
solutions we give in Fig. 3 the corresponding plot of p against the pressure P, which displays 
the pattern of a first order transition. The densities of the two stable phases are indicated 
by the light dashed lines in Fig. 2. Fig. 4 shows the T = 2°K and the T = Q°K isotherms 
in a standard P x p~^ diagram together with the appropriate Maxwell construction for 
T = 2°K. The condensed solution is the only one stable at T = 0°K. The pattern found for 
T — 2°K repeats itself for higher temperatures with no evidence of a critical temperature 
within the range allowed by the limitations of the effective theory. The stable isotherms for 
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temperatures of 0, 2, 3 and 4 "K are shown in Fig. 5. 

5 Concluding remarks 

In the gaussian approximation it is possible to renormahze a repulsive contact interaction 
using a procedure akin to the so called precarious renormalization scheme of Stevenson [pD|]. 
Despite being finite and exhibiting some desirable properties, such as an excitation spectrum 
of the phonon type, the condensate phase is thermodynamically unstable for all system 
densities. This is a manifestation of a corresponding instability poited out in a relativistic 
context by Bardeen and Moshe |T^. Solutions with no pairing {ak = 0), on the other hand. 



can be verified to correspond just to a free bose gas. We stress that thethermodynamic 
instability of the paired, condensed phase cannot be recognized just from the expression for 
the ground state energy, Eq. (27), but requires the use of the equation of state. 

As an alternative to such an unstable, renormalized theory we considered a cut-off de- 
pendent effective theory for a dilute system. Among the results of this theory we have a gap 
in the excitation spectrum in disagreement with the Hugenholz-Pines theorem[0, and a 
first order transition between a condensed phase and a non-condensed phase. The ocurence 
of the gap can be associated to the drastic reordering and infinite partial ressumation of the 
perturbative series implied in the Gaussian variational approximation. It should be noted 
in this connection that the variational calculation involving the truncated density Fq is de- 
signed to optimize the determination of the grand-potential, and that this does not imply the 
variational optimization of the thermal occupation probabilities from which the excitation 



spectrum is derived [^. On the other hand, the sensitivity of the excitation spectrum to 
the dynamical ingredients included in the calculation can be illustrated by the fact that, 
using the free quasi-boson truncation described in section 3 (which breaks the variational 
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self consistency of the calculation), the results obtained agree with the Hugenholtz-Pines 
theorem. 

We also remark that it is a well known tendency of mean field approximations to over- 
represent discontinuities in situations involving phase transitions, e.g. by predicting discon- 
tinuous behaviors in the case of finite systems. This recommends, as already observed in 



ref. p6[ on the basis of a reliability analysis of power-counting type, that caution should be 
exerted in the interpretation of the first order character obtained for the phase transition in 
this calculation. 

In conclusion, the introduction of quantum effects at the level of a paired mean field 
calculation for a non-ideal bose system gives results which differ even qualitatively from 
usual perturbative results. Comparison with experiment is still hampered by the fact that a 
very dilute though sufficiently non-ideal system has not yet been realized in practice. 
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Figure Captions 



FIG. 1. Cut-off dependence of the ground state energy per particle. See text for details. 

FIG. 2. Chemical potential as a funcion of number density for T=2°K and A = 3A ^ . See 
text for other parameter values. The highlighted values of fi and p correspond to the Maxwell 
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construction shown in Figs. 3 and 4. 

FIG. 3. Chemical potential as a function of pressure for the isotherm depicted in Fig. 2. 
The highlighted point indicates the first order phase transition. 

FIG. 4. Equation of state isotherms for T=0 and 2 "K showing the Maxwell construction 
identifying the phase coexistence region in the T—2 °K isotherm. 

FIG. 5. Equation of state stable isotherms for T=0, 2, 3 and 4 °K. 
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